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Abstract 

C^ , The stability of a Bose-Einstein condensed state of trapped ultra-cold 

atoms is investigated under the assumption of an attractive two-body 
and a repulsive three-body interaction. The Ginzburg-Pitaevskii-Gross 
' O ■ (GPG) nonlinear Schrodinger equation is extended to include an effective 

Ch ' potential dependent on the square of the density and solved numerically 

O , for the s— wave. The lowest frequency of the collective mode is deter- 

mined through the Fourier transform of the time dependent solution 
and its dependences on the number of atoms and the strength of the 
three-body force are studied. We show that the addition of three-body 

^ ' dynamics can allow the number of condensed atoms to increase con- 

OO \ siderably, even when the strength of the three-body force is very small 

compared with the strength of the two-body force. 



Q«^ ■ Recently, the theoretical research on Bose-Einstein condensation (BEC), a 

^\ , phenomenon predicted by Einstein more than 70 years ago, received consider- 

able support from the expedmental evidences of BEC in magnetically trapped 
2 ! weakly interacting atoms fl'D^B). The nature of the effective atom-atom in- 

teraction determines the stability of the condensed state: the two-body pseu- 
dopotential is repulsive for a positive s— wave atom-atom scattering length 
rH ■ and it is attractive for a negative scattering length p) . The ultra-cold trapped 

Q , atoms with repulsive two-body interaction undergoes a phase-transition to a 

O ' stable_Bose condensed state, irua number of cases found experimentally, as for 

^^Rb fJ), for ^^Na (p) and "^Li p). However, a condensed state of atoms with 
negative s— wave atom-atom scattering length would be unstable unless the 
. , number of atoms is small enough such that the stabilizing force provided by 

j^ ' the harmonic confinement in thatrap overcomes the attractive interaction, ^s 

found on theoretical grounds pQ). It was indeed observed in the ^Li gas p). 



"Contribution to the Internationl Workshop on Collective Excitations in Fermi and Base 
Systems, Serra Negra, Brazil, September 14-17, 1998. To be published by World Scientific, 
Singapore. 



for which the s— wave scattering length is a = (—14.5 ± 0.4) A, that the num- 
ber of allowed atoms in the Bose condensed state was limited to a maximum 
valueJDetween 650 and 1300, which is consistent with the mean-field predic- 
tion p) . An earlier experiment fl) suggested that the number of atoms in the 
condensate state was significantly larger than the theoretical predictions with 
two-body pseudopotential. This is consistent with an addition of a potential 
derived from three-body interaction, which can extend considerably the region 
of stability for the condensate even for a very small strength of the three-body 
force. 

f[panding our universe of possible effective interactions, it was reported in 
) that a sufficiently dilute and cold Bose gas exhibits similar three-body 
dynamics for both signs of the s— wave atom- atom scattering length. They 
concluded that the long-range three-body interaction between neutral atoms 
is effectively repulsive for either sign of the scattering length. Supposing that a 
repulsive three-body effective interaction is present in the atomic system, one 
can imagine that a stable condensate will be formed in the trap for a number 
of bosons large enough, such that thepthree-body interaction overcomes the 
attraction of the two-body interaction p) . 

In the present work, we investigate the competition between the attractive 
two-body interaction, originated from a negative two-atom s— wave scattering 
length, and a repulsive three-body effective interaction. (The latter can be 
originated by the existence of a weakly bound three-boson state, as it will be 
explained.) We show that, in a dilute gas, a small repulsive three-body force 
added to an attractive two-body interaction is able to stabilize the conden- 
sate beyond the critical number of atoms in the trap, found just with attrac- 
tive two-bndy force (EI). The Ginzburg-Pitaevskii-Gross nonlinear Schrodinger 
equation f3) is extended to include an effective potential dependent on the 
square of the density and solved numerically for s— wave. The stability is stud- 
ied using a weak perturbation probe and the time dependent equation is solved 
with the ground state as the starting point. The lowest frequency of the collec- 
tive mode is determined through the Fourier transform of the time dependent 
solution and its dependences on the number of atoms and the strength of the 
three-body force are studied. In this case, a signature for a repulsive three 
body force can appear with the possibility of two stable solutions, for a fixed 
number of atoms, even for a very small repulsive three body interaction. Also, 
in this case, one stable solution is found for a large number of atoms, beyond 
the maximum expected with only two-body interaction. 

In the following, we present the formalisHi, where the original Ginzburg- 
Pitaevskii-Gross (GPG) non- linear equation ^), which includes a term propor- 
tional to the density (two-body interaction) , is extended through the addition 



of a term proportional to the squared-density (three-body interaction). Next, 
after reducing such equation to dimensionless units, we study numerically the 
s— wave solution by varying the corresponding dimensionless parameters, which 
are related to the two-body scattering length, the strength of the three-body 
interaction and the number of atoms in the condensed state. As particularly 
observed in Ref. ^), to incorporate all two-body scattering processes in such 
many particle system, the two-body potential should be replaced by the many- 
body T—matrix. Usually, at very low energies, this is approximated by the 
two-body scattering matrix, which is directly proportional to the scattering 
length p). So, to obtain the desired equation, we first consider the effective 
Lagrangian, which describes the condensed wave-function in the Hartree ap- 
proxiniation, implying in the Ginzburg-Pitaevskii-Gross (GPG) energy func- 
tional ^): 
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In our description, the atomic trap is given by a rotationally symmetric har- 
monic potential, with angular frequency lo, and Ci gives the effective atom 
interactions up to three particles. 

The effective interaction Lagrangian for ultra-low temperature bosonic 
atoms, including two and three-body scattering at zero energy, is written as: 



£1 = -i /dVidV2dVidV2*t(r'i)*^(r'2)*(ri)*(r2) 
X (P12 |r(2)(0)| ri2) J3(^i + r-'a - n - r2) 
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where ri2 and R3 are the relative coordinates, given by fi2 — fi ~ -rz and 
R3=r3- (ri + r2)/2; and i?i23 = (n + r 2 + fg). T(3)(0) and T^^\o) are 
the corresponding three-body T—matrix and two-body T—matrix for the pair 
jk, which are evaluated at zero-energy. The two-body T—matrix for each pair 
(jk) is subtracted from T^^\0) to avoid double counting and Ki is the kinetic 
energy operator for particle i. 



We can approximate the above effective interaction Lagrangian at low den- 
sities by averaging the T— matrices over the relative coordinates, considering 
that the thermal wave-length is much greater than the characteristic interac- 
tion distances. 
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The integrations of the T-matrices over the relative coordinates gives the zero 
momentum matrix elements: 



dM2dVi2(r^2|r(2)(0)|ri2) = (27r)3^pi2 = 0|r(2)(0) 



P12 = 



Airfi^a 



(4) 



where a is the two-body scattering length. For the connected three-body 
T— matrix, we have 
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where A3 is the strength of the three-body effective interaction. 

The nonlinear Schrodinger equation, which describes the condensed wave- 
function in the mean-field approximation, is obtained from the effective La- 
grangian given in Eq. (nf). By considering the inter|action in Eq. (H) and the 
definitions in Eqs. (0) and (0), it can be written as f^) 



2m 2 



|*(f)p + A3|*(r)|^ 



*(f) = M*(f) , (6) 



where fj, is the chemical potential, fixed by the number N of atoms in the 
condensed state: 

'"dV|*(f)p ^ N . (7) 

The physical scales present in Eq. (0) can be easily recognized by working with 
dimensionless equations. By rescaling Eq.(p) for the s— wave solution, we ob- 
tain 



dx^ 4 x^ x^ 



^x) = P^x) , (8) 



where x = sj2vrujjj'h r and $(x) = •\/87r|a| r'^{r). The dimensionless pa- 
rameters, related to the chemical potential and the three-body strength are, 
respectively, given by 



/3 = - — and 53 = X^fiuj 
Tioj 



Anh'^a 



(9) 



and the normalization for ^{x), obtained from Eq. (M), defines a number 
related to the number of atoms N: 



/■OO fly 

/ dxmx)f =n, where n = 2N\a\J-^. (10) 

The boundary conditions in Eq.(P) are given by p) 

$(0) = and $(00) - Cexp[-xV4-l-(/3-l/2)ln(x)]. (11) 

The numerical solutions of Eq. (0) are obtained for several values of /3, 
using three values of 173 to characterize the solutions. We have used the Runge- 
Kutta (RK) and "shooting" method to obtain the corresponding solutions in 
each case. The stability assignment was made by studying the corresponding 
time dependent Scrodinger equation, using the Crank-Nicolson (CN) method 
[see Ref. fj)]. The numerical procedure to determine such stability was done 
in the following way: when applying the CN method, we started by using 
the static solution obtained from the RK method and observed if the modu- 
lus of the wave function remained constant. If this was occuring for a long 
period of time (of about 500 units of dimesionless time r — ujt) the solution 
was considered stable, otherwise unstable. In another procedure we added a 
weak perturbation to the potencial. We examined the time evolution of a se- 
lected point of the wave- function. The lowest collective ospilations (uJcoi) were 
determined by using the Fourier transform of such result ^J) . 

Figures 1-4 presents our main results. 




Fig. 1 The chemical potential (/3) in units of Tiuj as a function of n for three 
values of the non-dimensional strength of the repulsive three-body effective interac- 
tion (gs). fl3 = (dashed line); gs, = 0.016 (solid line); gs = 0.03 (dot-dashed line) 
[see Eqs. (bl) and (|lO|)] . In the inset, for gs = 0.016, the solutions linking A and B 
are unstable. 
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Fig. 2 Frequencies of the lowest collective excitations (with I — 0) as a, function 
of n for three values of the dimensionless strength gs The points A and B correspond 
to the points shown in the inset of Fig.l, where the frequencies go to zero. 
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Fig. 3 Frequencies of the lowest collective excitations (/ = 0), as a function of 
/3 for three values of the nondimensional strength gg. The points A and B correspond 
to the points shown in the inset of Fig.l, where the frequencies go to zero. 




Fig. 4 Total energy of the Bose condensate, given in units of hco/ I 2\a\^ {2rauj jTi 
is shown as a function of n, for three values of the nondimensional strength gs. The 
line linking A and B in the inset correspond to the unstable region shown in Fig. 1. 



In Figure 1 we present (3 as function of n for 173 equal to 0, 0.016 and 
0.03. Our calculation for 53 = agrees with the result presented in Ref. p), 
with the maximum number of atoms limited to rimax ~ 1.62 |j. As we can 
see, for n < Umax two solutions are possible, the lower being unstable. For 
(73 slightly higher than a new pattern appears, as one can see for 173 = 0.016 
presented in the figure. We can divide the figure, for 173 = 0.016 (for example), 
in three regions; the first, stable, goes from f3 = 1.5 until point A (presented in 
the inset); an unstable region appears from the point A till B; another stable 
region goes from B towards n — > 00. 

Figure 2 shows the collective frequencies lUcoI as function of n for the first 
mode {I = 0). The^ solutions corresponding to 53 = agree well with the 
ones given in Ref. Ipl); and they loose stability as ujcoi -^ 0. By this criteria, 
we determined the regions of stability for 53 = 0.016. For 53 = 0.03 all the 
solutions are stable. The collective frequencies are also shown as a function 
of /3 in Figure 3. The total energy is shown in Figure 4. jTjhe case (73 = 
reproduces the results already found by Houbiers and Stoof pJ) , including the 
unstable (higher) solutions. For 53 = 0.016 the part of the plot linking points A 
and B is unstable, corresponding to the unstable region of Figure 1. Otherwise 
it is stable. Finally, for 53 = 0.03, the function of the energy in terms of n is 
always single valued and stable. 

The relevance of the three-body interaction was neveLcmphasized in BEC 
studies. However, in 1985, it was pointed out in Ref. Ipl) that an easier ex- 
perimental approach to probe density fluctuations is to consider an observable 
directly sensitive to the probability of finding three atoms near each other, 
which will correspond to the loss rate of atoms due to three-body recombina- 
tion. Suck, a t|hree-bodii-|recombination rate in BEC, was considered recently 
in Refs. (ta), P) and ^) It was shown in Ref. pl) that the three-body re- 
combination coefficient of ultcacold atoms to a weakly bound s level goes to 
infinity in the Efimov limit ^). The Efimov limit is a particularly interesting 
three-body effect, which happens when the two-body scattering length is very 
large (positive or negative). In this case, with the two boson energy close to 
zero, the three-boson system presents an increasing number of loosely bound 
three body states, which haveuJarge spatial extension and do not depend on 
the details of the interaction P3). In such limit, when the three-body bind- 
ing is close to zero, the strength A3 can be large enough to give a sensible 
value for 53. So_-Our main motivation here was to provide an extension to the 
GPG equation ^), which considers a three-body interaction and, in this way, 
provides the framework for a numerical investigation of the relevance of three- 
body interaction in Bose-Einstcin condensation. In this letter, we presented 

'' Our n is equal to \Ci^f\ of Ref. |l). 



results for the chemieal potential, fi, the number of atoms in the eondensed 
state (N), and frequency of the lowest collective mode for a range of values 
of A3, expressed respectively in terms of the dimensionless parameters /?, n, 
^coil^trap and 53. The determination of A3 has to be considered in particular 
three-boson systems. 

To summarize, our calculation presents, at the mean-field level, the con- 
sequences of a repulsive three-body effective potential added to an attractive 
two-body interaction for the Bose condensed atomic state. We propose a gen- 
eral framework for discussing such phenomenon. As we have shown with the 
present approach, a slightly repulsive three body interaction provides stable 
solutions for higher number of atoms where their existence is forbidden when 
considering only two body potential. This can be a signature of the presence 
of an effective repulsive three body interaction in Bose condensates. 
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